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ABSTRACT 



In order to investigate the origin of the excess of strong Mgn systems towards GRB afterglows as compared to QSO sightlines, 
we have measured the incidence of Mgn absorbers towards a third class of objects: the Blazars. This class includes the BL Lac 
object population for which a tentative excess of Mgn systems had already been reported. We observed with FORS1 at the ESO- 



VLT 42 Blazars with an emission redshift 0.8 < 



1.9, to which we added the three high z northern objects belonging to the 



Uy BL Lac sample. We detect 32 Mgn absorbers in the redshift range 0.35-1.45, leading to an excess in the incidence of Mgn 
absorbers compared to that measured towards QSOs by a factor ~ 2, detected at 3cr. The amplitude of the effect is similar to that 
found along GRB sightlines. Our analysis provides a new piece of evidence that the observed incidence of Mgn absorbers might 
depend on the type of background source. In front of Blazars, the excess is apparent for both 'strong' (vv,(2796)> 1.0 A) and weaker 
(0.3 < w r (2796)< 1.0 A) Mgn systems. The dependence on velocity separation with respect to the background Blazars indicates, at 
the ~ 1.5cr level, a potential excess for/? = v/c ~ 0.1. We show that biases involving dust extinction or gravitational amplification 
are not likely to notably affect the incidence of Mg n systems towards Blazars. Finally we discuss the physical conditions required 
for these absorbers to be gas entrained by the powerful Blazar jets. More realistic numerical modelling of jet-ambient gas interaction 
is required to reach any firm conclusions as well as repeat observations at high spectral resolution of strong Mg n absorbers towards 
Blazars in both high and low states. 

Key words. BL Lacertae objects: general - quasars: absorption lines 



1. Introduction 

In the 1980's, the observation of large QSO samples allowed 
Young et al. (11982b and a decade later Steidel & Sargent d 19921 1 
to establish that C iv and Mg n absorption systems arise in inter- 
vening material, randomly distributed along the QSO sightlines. 
This conclusion was based on two main arguments: i) the num- 
ber of systems detected in individual sightlines is Poissonian; ii) 
the distribution of systems with respect to comoving coordinates 
is uniform. Their Mg n sample did not reveal any cosmological 
evolution. 

To firmly establish the nature of the lower redshift Mg n ab- 
sorption systems, searches for galaxies associated with these ab- 
sorptions have been conducted since the late 1980's. For Mgn 
systems at z < 1, results of these identification surveys revealed 
that normal field galaxies with luminosities L > 0.1L* are all 
potentially Mg n absorbers, and the mean extent of their gaseous 
halos is ~65 kpc, assuming Hq = 70 km s _ Mpc~' (Bergeron 
& Boisse [HHJ Steidel et al. [19941 Steidel [19951 The conclu- 
sion of these searches was that in studying Mgn systems, we 
get information on the cosmic history of normal field galaxies 
(metallicity, dust content, dynamics, star formation). 

Our current understanding of absorption line systems thus 
implicitly assumes that their properties do not depend on the 
nature of the background sources used as probes (except at 



* Based on observations collected at the European Southern 
Observatory, ESO, the VLT/Kueyen telescope, Paranal, Chile, in the 
framework of programs 080.A-0276 and 081 A-0193. 



Zabs - Zem)- However, Prochter et al. (2006a) recently reported an 
excess of strong (rest equivalent width w r (Mg n2796) > 1.0 A) 
Mg ii absorbers towards gamma-ray bursts (GRBs) with respect 
to the incidence measured in QSO spectra. These authors found 
that nearly every GRB sightline exhibits at least one absorber 
((z) = 1.1) whereas the incidence of strong Mgn QSO ab- 
sorbers is about four times smaller. This excess is confirmed by 
Sudilovsky et al. 12007) for a small, homogeneous sample of 
GRB VLT-UVES high-resolution spectra, and more recently by 
Tejos et al. ( 120091 1 and Vergani et al. (|2009). Based on their im- 
proved statistics, these latter authors find the excess to be smaller 
than initially reported, by a factor of 2.1 + 0.6 instead of ~4. 
There is no such excess for the higher redshift C iv absorbers 
(small overlapping samples of Sudilovsky et al. 2007 and Tejos 
et al. 120071 1. which trace more diffuse and hotter gas than Mg n, 
have larger cross-sections and, therefore, largely probe a differ- 
ent population than strong Mg n absorbers. 

Studies of systems towards GRBs naturally focused on the 
strongest systems, the only ones that could be easily detected and 
thus for which homogeneous results could be obtained (since 
they fade away rapidly, GRB afterglows are generally faint when 
observed at good spectral resolution). Initially, the incidence of 
strong Mg n systems in QSO spectra was poorly known because 
their number density rapidly decreases with increasing w, (Mg ii) 
(e.g. Steidel & Sargent [1992l . and only the SDSS data releases 
provided a robust measurement to be used for the purpose of 
comparison with the GRB results. As time elapsed, the num- 
ber of afterglows observed at a resolution and signal-to-noise ra- 
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tio (S/N) sufficient to detect weaker systems increased, allowing 
Vergani et al. (20091 to measure the incidence of Mg n systems 
with 0.3 < w r (Mgn2796)< 1.0 A and to show that the latter do 
not display an excess towards GRBs as compared to QSOs. Thus 
the GRB excess is only significant for strong Mgn systems. 

Several physical effects have been proposed to explain the 
GRB-QSO absorber discrepancy: 1) obscuration by dust of 
background QSOs, 2) intrinsic nature of the GRB absorbing gas, 
3) gravitational lensing of GRBs by Mg n absorbers, 4) absorbers 
with small cores of a similar size as GRB and QSO beam sizes 
(Prochter et al. l2006al Frank et a l.l2007l Hao et al. [2007l . 

In their study, Prochter et al. d2006al) mention a possible ex- 
cess of the same kind for strong Mg n systems towards BL Lac 
objects, as inferred by Stocke & Rector ( 119971 ). BL Lac objects 
can be of great interest in our context because these targets are 
mainly radio-selected and thus should not suffer from any ex- 
tinction bias. Furthermore, the timescale and amplitude of their 
optical variations is intermediate as compared to those of QSOs 
and GRBs. Finally, relativistic boosting is thought to be present 
in BL Lac objects, as in GRBs. 

Since the result obtained by Stocke & Rector was based on 
relatively few BL Lac spectra, not all with good S/N, and thus 
of poor statistical significance, we have undertaken an observa- 
tional programme devoted to the measurement of dN/dz towards 
Blazars for strong and weak Mg n systems. Indeed BL Lacs be- 
long to the broader Blazar class characterized by violent variabil- 
ity, a compact flat-spectrum radio source, and a smooth optical- 
near infrared continuum. Since the available data on absorption 
systems in BL Lac objects with the appropriate resolution and 
S/N has been limited so far, our survey can also provide valu- 
able results on systems with an intrinsic nature, to be compared 
with the detailed information provided by the SDSS data on as- 
sociated systems in QSO spectra. One might expect for instance 
differences related to the nature of the host galaxies or to the 
presence of interstellar gas entrained by the jet. 

The paper is organised as follows. In Sect. 2, we describe 
how targets were selected. We present the observations and raw 
results in Sect. 3. Results on the incidence of intervening and as- 
sociated Mg ii systems are given in Sect. 4. In Sect. 5, we com- 
pare our results with those obtained for GRBs and QSOs and 
investigate possible interpretations of the excess of Mg n Blazar 
absorption systems. We give our summary in Sect. 6. 



2. The sample 

We aim at measuring dN/dz for Mgn intervening absorption 
systems towards BL Lac objects with both strong and weak 
rest equivalent width, in order to check whether the excess, if 
real, depends on the strength of the Mgn absorption doublet. 
This is not an easy task for several reasons. First, our project 
requires targets with known emission redshifts to properly es- 
timate the redshift interval probed. Second, BL Lacs and more 
generally Blazars exhibit large, unpredictible magnitude fluctu- 
ations, making it difficult to obtain a homogeneous data set. As 
in previous studies, we shall consider for our statistical analyses 
that systems with a relative velocity Av > 5000 km s _I from the 
source are of intervening nature while the others are considered 
as associated. 



2.1. Target selection 

Since, at low redshift (z < 1.2), there is a strong cosmic 
evolution of the number density of QSO Mg n absorbers with 



Table 1. The sample. 



target 


class 


coordinates 




grism 


At/n 


common name 




2000 






min 


PKS 0057-338 


BL 


/ \ 1 t \/ \ t \t \ 1 \ T)OT1A 

010009.0-333730 


0.875 


B600 


120/3 


PKS 0208-512 ' 


BL 


021046.2-510101 


1.003 


B600 


80/2 


PKS 0215+015 7 


BL 


021748. 9+014449 


1.715 


V600 


60/2 


2QZ 


opt 


023405.5-301519 


1.690 


V600 


80/2 


PKS 0235+164 ' 


BL 


023838.9+163659 


0.940 


B600 


60/3 


SDSS 


opt 


f\r\ A 1 £ S A , f\f\ Ai-\ e 1 

024156.4+004351 


0.989 


B600 


120/3 


FKS 0256+0 /5 


BL 


U2592/.0+U/4/39 


0.893 


B600 


80/2 


PKS 0332-403 ' 


BL 


033413.6-400825 


1.351™ 


V600 


115/3 


PKS 0420-014 ' 


FS 


042315.8-012033 


0.915 


B600 


80/3 


r> \/ ci CiA^tC TO f\RS V 

FKS 042o-3o0^ y 


BL 


f\A^O Af\ A QTCtm 

U42840.4-3 /5ol9 


1.110 


B600 


80/2 


PKS 0454-234 ' 


BL 


045703.1-232452 


1.003 


B600 


90/3 


n T r n AT'IH AAA RS "vC 

PKS 0537-441 yL 


BL 


053850.3-440508 


0.890 


B600 


20/4 


r»T.?"o m a i r\f\/' 

PKS 0743-006 


FS 


074554.0-004417 


0.994 


B600 


100/3 


PKS 0808+019 


BL 


081126.7+014652 


1.145 


B600 


80/2 


PKS 0823-223 />Gr 


BL 


082601.5-223027 


>0.910 a 


B600 


20/2 


PKS 0906+015 


FS 


f\f\f\r\ 1 r\ 1 . a i o i i r 

090910.1+012135 


1.022 


B600 


25/1 


SDSS 


BL 


f\C\ A /~\ C ^ O f\f\ A TAT 

094257.8-004705 


1.362 


V600 


80/2 










B600 


40/1 


SDSS 


BL 


094827.0+083940 


1.489 


V600 


80/2 


SDSS 


BL 


100959.6+014533 


> 1.085" 


V600 


120/3 


PKS 1008+013 


BL 


1 A 1 1 1 C z' . f\ 1 CX/' A o 

101 1 15.6+010642 


1.275™ 


V600 


80/2 


PKS 1130+008 


BL 


1 1 O ^ A C /" . AA1 /in 

113245.6+003427 


1.234 


V600 


60/2 


PKS 1 144-3 /9 ' 


BL 


1 14/U1.3— Jol21 1 


1.049 


B600 


40/2 


SDSS 


BL 


124533.7+022825 


> 1.096" 


B600 


80/2 










V600 


80/2 


PKS 1250-330 


BL 


125258.3-331959 


0.856 


B600 


120/3 


FIRST" 


BL 


133859.0+1 15316 


>1.589 b 


V600 


80/2 


PKS 1406-076 yC 


FS 


140856.4-075226 


1.500 


V600 


80/2 


PKS 1407+022 7,0 


BL 


141004.6+020306 


1.253™ 


V600 


80/2 


SDSS 


BL 


141927.4+044513 


> 1.684" 


V600 


80/2 










R600 


30/1 


FKS 1424-415 ' 


FS 


142/5o.2-42Uol9 


1.522 


V600 


57/2 


PKS 1519-273 


BL 


152237.6-273010 


1.294 


V600 


115/3 


SUSS 


BL 


o4M5. /+U0323j 


1.031 


V60U 


OU/2 


FKS 1 /41-035 ' 


FS 


1 /435o.o-U350U4 


1.054 


B600 


80/2 


PKS 1953-325 


FS 


195659.4-322546 


1.242 


V600 


120/3 


PKS 2012-017 re 


BL 


201515.1-013732 


>0.940 c 


B600 


80/2 


PKS 2029+ 121^ 


BL 


203155.0+121941 


1.215 


V600 


80/2 


PKS 2131-021 MC 


BL 


213410.3-015317 


1.285 


V600 


80/2 


MH2133-449 />G 


BL 


213618.3-444348 


>0.980 c 


B600 


120/3 


EQ 2207.3+004 7,0 


BL 


220719.7+004157 


1.892 


V600 


80/2 


2QZ 


opt 


221450.1-293225 


1.636 


V600 


80/2 


PKS 2223-114 PG 


BL 


222543.7-111341 


0.997 


V600 


80/2 


PKS 2223-052° 


BL 


222547.2-045701 


1.404 


V600 


80/2 


PKS 2308-109 


BL 


231116.9-103849 


1.529 


V600 


80/2 


B2 0218+35'"''' 


BL 


022105.4+355615 


0.944 


N/A 


N/A 


S5 0454+844* 5 


BL 


050842.4+843204 


> 1.340 


N/A 


N/A 


B2 1308+32^^ 


BL 


131028.7+322044 


0.997 


N/A 


N/A 



class - BL: BL Lac object, opt: optical BL Lac candidate, FS: Flat 

Spectrum Radio QSO (FSRQ). 
n Number of exposures. 

RS In the 1 Jy sample of Rector & Stocke (2001). 
PG In the BL Lac sample of Padovani & Giommi (1995) and not in RS. 
pl In the BL Lac list kindly provided by R. Plotkin prior to publication. 
7 Gamma-ray emitter. 
c Superluminal motions. 
N New emission redshift. 
" Upper limit derived from Mg n absorption. 
b Upper limit derived from C iv absorption. 
c Upper limit derived from lack of Ca n absorption from the host 
galaxy. 
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w r (Mgn2796) > 1 A (Nestor et al. l20531 Prochter et al. l2006bt . 
we selected targets with an emission redshift z em > 0.8 to opti- 
mize the efficiency of our spectroscopic observations. This also 
yields a redshift range comparable to that investigated for GRBs 
by Prochter et al. (2006 a|). We restricted further the redshift 
range to z em < 1.8 due to the paucity and faintness of Blazars 
at higher redshift. 

Our sample comprises the seven 1 Jy BL Lacs of Rector & 
Stocke (T200T1) at z em > 0.8 and DEC < +20deg to: 1) confirm 
the strength of the Mg n absorption lines for spectra of low S/N, 
2) detect weaker Mg n, 3) complement the available information 
on emission redshifts. The same selection criteria were applied 
to the BL Lac sample of Padovani & Giommi i 19951 1 (nine addi- 
tional targets). We have taken into account recent emission red- 
shift measurements, including upper limits (Sbarufatti et al. 2005 
& 120061) . The latter are derived from either the Mg n absorption 
system of highest redshift or the absence of the Ca n absorption 
doublet from the host galaxy in the BL Lac spectrum. We have 
re-observed almost all the 0.9 < z e m < 1-8 BL Lacs studied 
by these authors, since the Mg n absorption doublets are not re- 
solved in their low resolution FORS1 data (FWHM = 15-20 A 
with the 300V+I grism). 

To increase our sample, we searched for targets, at 0.8 < 
Zem < L8 and DEC < +20deg, of the Blazar class in the 
SIMBAD and NED databases. In addition to BL Lacs, the Blazar 
class also includes Flat Spectrum Radio QSOs (FSRQs) which 
exhibit broad emission signatures (Angel & Stockman 119801 1. 
Blazars are thought to have jets pointing towards us with asso- 
ciated emission beamed by relativistic motions. When in high 
state, the optical emission can be highly polarized (e.g. Impey 
& Tapia 1 1 9881 & 1 1990b . When in low state, the isotropic ther- 
mal emission from the active nucleus and its associated broad 
emission line region can be detectable or even dominate over the 
Doppler boosted emission from the jet. The magnitudes given 
in SIMBAD only reflect the state of the object at a given epoch; 
nevertheless, we introduced a V magnitude cut of 2 1 .0 to exclude 
(not always successfully) very faint sources. 

This SIMBAD and NED searches provided most of the re- 
maining targets of our sample. It includes optically selected BL 
Lacs among t he 2DF AATQSO redshift s urvey a nd the SDSS 
(Londish et al. [20021 & [20071 Collinge et al. 120531 . Some of the 
SDSS BL Lacs have high optical polarization which confirms the 
synchroton nature of the optical emission (Smith et al. l2007l >. 

Finally, we considered information kindly provided by 
Plotkin on seven 'southern' high-redshift BL Lacs prior to pub- 
lication (Plotkin et al. 12008b of their BL Lac sample selected 
from a combination of the SDSS Data Release 5 (DR5) and the 
VLA FIRST survey. Only two of them, at z em < L8, were not 
already in our SIMBAD/NED list (see Table G}; both were ob- 
served. The upper limit on z em for one of them is derived from a 
Civ absorption doublet detected in the SDSS spectrum. 

The resulting sample comprises 5 1 sources for our observa- 
tions with FORS1 scheduled in November 2007 and May 2008. 
Due to crowding in right ascension (RA) around 11- 12h, we 
could only observe 42 of them; these were chosen according to 
their RA and DEC only, disregarding any previously available 
information on absorption systems. To consistently take into ac- 
count all the 1 Jy sources of Rector & Stocke ( 120011 1 at z em > 0.8, 
we included the three northern sources which have published, 
reliable spectroscopic information on the presence or lack of 
strong Mg n systems (Miller et al. 119781 Stocke & Rector 1 19971 
Cohen etal. 120031 . 

Thus, there are 45 sources at 0.8 < z em < 1.8 in our final 
sample of which ten BL Lacs from the 1 Jy sources of Rector 



& Stocke (2001 ). Our sample comprises 35 BL Lacs, 7 FSRQs 
and 3 optical BL Lac candidates. Among the radio-loud Blazars, 
there are 15 gamma-ray sources, this emission being an impor- 
tant component of the relaticistic jet (Mukherjee et al. 119971 
Abdo et al. |2009l >. Superluminal motions are detected in 1 1 ra- 
dio Blazars (Piner et al. 120061 120071 Fan et al. |2009l and, of 
those, eight are indeed also known gamma-ray emitters. Our 42 
observed sources are presented in Table[T]as well as information 
on the three northern BL Lacs. We now discuss possible biases 
affecting this sample. 

2.2. Potential biases 

In order to obtain a reliable measurement of dN/dz, one must 
carefully identify potential selection effects that might bias the 
estimate towards larger (or smaller) values. 

One first difficulty is related to the fact that for some sources, 
no emission lines had been detected or identified at the time of 
our sample selection and only a lower limit on z em was provided 
by a Mg n system. Including such systems in the statistics would 
clearly bias dN/dz towards a too high value since they might in 
fact be associated instead of intervening. To avoid this bias, we 
discard these possibly associated Mg n systems, if our observa- 
tions do not reveal new emission lines providing a firm source 
redshift. For these cases, we include in the analysis only the 
lower redshift part of the sightline (i.e. beyond 5000 km s _1 of 
the absorption redshift). This is a conservative approach, since 
some of the 'possibly associated' Mg n systems could be inter- 
vening. Four of our targets have z em constrained by a strong Mg n 
absorber, and one by a high redshift C iv absorber without asso- 
ciated Mg ii. On the opposite, if our spectrum provides a new 
firm redshift for the source, all intervening Mg n systems are in- 
cluded. 

A second concern involves sources for which intervening 
systems were already known to be present. Ideally, one should 
observe for our programme only sources for which no infor- 
mation at all was available on the presence or absence of ab- 
sorption systems. Unfortunately, this would severely restrict the 
number of targets, which is already small. To circumvent this 
difficulty, we performed our selection exclusively on the basis of 
emission properties. If those targets with already known systems 
were excluded, this would bias dN/dz towards lower values. 
Although spectra exist for sources like PKS 0235+164 and PKS 
0215+015, we nevertheless reobserved them to obtain either a 
more complete wavelength coverage or a uniform w r (Mg n2796) 
limit. Further, this offers the opportunity to investigate temporal 
variations in these systems over a reasonably large time interval, 
since some of them were initially detected at high S/N about 25 
years ago (for PKS 0215+015 see Bergeron & D'Odorico [T986l 
we will present our results on possible variations of Mg n and 
Fe ii absorptions in a forthcoming paper). 

Finally, given the broad property range of the Blazar class, 
it remains possible that some targets be in fact QSOs. Since our 
main purpose is to test whether the number density of strong 
Mg ii systems is larger towards Blazars, this is however not a real 
problem because contamination of our sample by QSOs would 
then lower our dN/dz estimate. 

3. Observations 

3.1 . Observations and data reduction 

The observations were performed with the FOcal Reducer and 
low-dispersion Spectrograph (FORS1) at the ESO VLT UT2 
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Table 2. Mgn absorption systems. 



target 




emission lines 




vy r (2796,2803)" 
A A 


w r (2586,2600)" 
A A 


associated absorptions 


0100-3337 


0.875 


Cm],Mg ii 


0.6810 


0.36,0.36 


nd,nd 


- 


0217+0144 


1.715 


Cm] 


0.7557 


0.12,0.06 


nd,nd 


- 








0.9685 


0.22,0.15 


nd,nd 


- 








1.3439 


1.92,1.83 


0.75,1.39 


MgI,Fen, Mnn,Znn,Crn 


0234-3015 


1.690 


Aim], Cm] 


1.2361 


0.11,0.09 


nd,nd 


- 


0238+1636 


0.940 


Cm],Mg ii 


0.5245 5R 


2.26,2.12 


0.89,1.63 


MgI,Fen,Mnn 








0.8522 s * 


0.43,0.23 


nd,0.09 


Fen 


0241+0043 


0.989 


Cm],Mg ii 


0.6310 


0.27,0.18 


nd,nd 


- 








0.7753 


1.06,0.81 


0.21,0.39 


Fen 








0.9790 A 


0.64,0.38 


nd,nd 


Aim 


0334-4008 


1.351" 


Mg ii 


1.0791 


0.20,0.13 


nd,nd 


- 








1.2083 


0.84,0.77 


0.16,0.31 


Fen 


0423-0120 


0.915 


Cm],Mg ii 


0.6338 


0.88,0.70 


0.05,0.14 


Mgl, Fen 


0428-3756 


1.110 


Mgn 


0.5592 


0.93,0.72 


0.19,0.42 


MgI,Fen 








1.0283 s * 


0.56,0.41 


nd,0.08 


Fe ii, Al m 


0457-2324 


1.003 


Mg ii 


0.8922 


2.20,1.90 


1.38,1,87 


MgI,Fe ii, Mn n,Zn n,Crn,Al m 


0538-4405 


0.890 


Mgn 


0.6725 


0.065,0.041" 


nd,nd 


- 


0745-0044 


0.994 


Cm],Mg ii 


0.7979 


0.095,0.056 


nd,nd 


- 


0826-2230 


>0.910 


- 


0.7057 


0.11,0.06 


nd,nd 


- 








0.9107 A? 


1.28,0.92 


0.19,0.41 


MgI,Fen 


0909+0121 


1.022 


Cm], Mgn 


0.5369 


0.41,0.31 


nd,0.14 


Fen 


0942-0047 


1.362 


Mg n,broad abs 


0.8182 


1.58,1.28 


0.70,0.96 


MgI,Fe n,Mn n,Zn n,Cr n 








1.0231 


0.38,0.22 


nd,nd 


- 


0948+0839 


1.489 


Cm],Mg ii 


1.0763 


0.76,0.63 


nd,0.26 


Fen 








1.3273 


4.07,3.72 


1.80,2.80 


MgI,Fen,Mnn 








1.4247 


1.08,0.62 


0.16,0.37: 


Fen 


1009+0145 


> 1.085 


- 


1.085 1 A? 


1.23,0.90 


0.18,0.32 


Fen 


1147-3812 


1.049 


Cm],Cii],Mgn 


0.3750 


0.27,0.20 


- 


- 


1245+0228 


> 1.096 


- 


1.0116 


0.60,0.48 


0.26,0.44 


Fen 








1.0946 A? 


2.30,2.27 


1.82,2.04 


MgI,Fe n,Mn n,Zn n, Cr n,Al n,Al m,Si n 


1408-0752 


1.500 


Cm],Cii],Mgn 


1.2753 


2.08,1.90 


1.17,1.53 


MgI,Fe n,Mn n,Zn n,Cr ii 








1.2913 


0.27,0.26 


nd,0.06: 


Fen 


1410+0203 


1.253" 


Cn],Mgn 


1.1123 


0.79,0.56 


nd,0.08: 


Fen 


1419+0445 


> 1.684 


- 


0.8987 


0.73,0.91* 


0.21,0.35 


Fen 








1.0849 


0.47,0.37 


0.13, b 


Fen 








1.1039 


0.52,0.25 


nd,nd 


- 








1.1834 


0.12,0.09 


nd,nd 


- 








1.2722 


0.62,0.58 


b ,0.49 


Fen 








1.6832 A? 


2.60,1.89 


0.49,0.80 


Fen 


1427-4206 


1.522 


Cm],Cii],Mgn 


1.0432 


0.23,0.12: 


nd,nd 


- 








1.0907 


0.62,0.37 


nd,0.09 


Fen 


1522-2730 


1.294 


Mgn 


1.2847 A 


0.25,0.22 


nd,0.06 


Fen 


1743-0350 


1.054 


Cm],Cn],Mgn 


0.2527 


0.64,0.49: 


- 










0.5293 


1.05,0.80 


0.37,0.65* 


MgI,Fen 








0.9077 


2.46,2.05 


0.89,1.32 


MgI,Fe n,Mnn 


IVjo— iZlj 


1 .Z4Z 


Mg ii 


U.OZ^O 


n n on 












1.0660 


3.53,3.49 


1.77,2.70 


Mgl'\Fen,Mnn 


2031+1219 


1.215 


Mg ii 


1.1158 s * 


1.29,1.16 


0.71,1.12 


MgI,Fen,Mnn 


2134-0153 


1.285 


Mg ii 


1.2458 A 


0.22,0.15 


nd,nd 




2136-4443 


>0.980 




0.5211 


0.51,0.39 


0.12:,0.26 


MgI,Fen 


2214-2932 


1.636 


Alm,Cm] 


1.5070 


0.38,0.33 


nd,0.17 


Fen 


2225-0457 


1.404 


Cm],Cn],[Neiv],Mgn 


0.8458 


0.54,0.43 


0.05,0.18 


Fen 








0.9782 


0.10,0.07 


nd,nd 




2311-1038 


1.529 


Cm],Cn] 


1.5242 A 


3.61,3.09 


1.37,2.34 


MgI,Feii,Alm 


0221+3556 


0.944 


Mg ii 


0.6850 s * 


1.78,1.72 


na,na 


Fe n,MgI 



" Representative uncertainties are given in Sect. 3.2. 
N New emission redshift. 



Strong Mgn system in Stocke & Rector (1997). 
A Associated system. 
A? Possibly associated system. 

" Target in very bright state. The optical flux is known to vary by a factor of at least 35. 
b Blended with an other absorption, 
nd : not detected, 
na : not available. 
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Fig. 1. Redshift path density for the intervening Mg n systems 
towards Blazars: strong (w r (2796) > 1.0 A: solid black curve) 
and weak ones (0.3 < w, (2796) < 1.0 A: thinner red curve). 



(Kueyen) telescope in November 2007 and May 2008 (Visitor 
Mode). During these observing runs, the seeing was in the range 
0.5" to 1.3", with an average of 0.8". In the long-slit mode, 
the selected B600 and V600 grisms covered the 3405-6095 A 
and 4435-7390 A spectral ranges, respectively, and yielded a 
spectral resolution, measured on our target spectra, of FWHM= 
5.0 + 0.5 A (or 3.4 + 0.34 binned pixels) for our 1.0" slit width. 
Details on the exposure times and grisms are given in Table Q] 
Data reduction was performed using MIDAS (08FEBpl 1.1). The 
double exposures for each target enabled an adequate correc- 
tion of cosmic ray contamination. Standard spectroscopic stars 
(HILT 600, LTT 9239 and for both observing runs LTT 3218) 
were observed with a 5.0" slit width. The seeing although gen- 
eral good was nevertheless fairly variable and thus only a relative 
flux calibration could be performed. 

3.2. Redshift and equivalent width estimates 

Among the 42 Blazars observed with FORS1, we detect Mgn 
and/or C m] or C n] emission lines in 35 sources. All the detected 
emission lines are given in Table [2] Some of our z em estimates 
differ from those those already published, including SDSS ones, 
as our spectra are of higher S/N. We could determine new emis- 
sion redshift (labelled N in Tables Q] and [2]) for three Blazars, of 
which one without any Mgn absorption system. Their spectra 
are presented in Appendix A. 

Mgn absorption doublets are detected in 31 Blazars ob- 
served with FORS1 plus in one of the three northern 1 Jy BL 
Lac. Their absorption redshift, z a b s , and rest equivalent widths 
are listed in Table [2] We also give the rest equivalent widths 
of their associated Fen2586,2600 doublets together with a list 
of all the other detected absorptions. Uncertainties on w b S val- 
ues vary depending on the spectral region and target magnitude. 
Representative values are in the range 0.04 - 0. 12 A (correspond- 
ing to relative uncertainties of about 2 to 10% depending on line 
strength) for intermediate S/N values of 50 to 100. In our es- 
timate, we include both pixel-to-pixel noise and uncertainty in 



the continuum placement; the latter is dominant when the S/N is 
larger than about 100. In some cases, the Blazars were in a very 
high state with an optical flux fully dominated by synchrotron 
emission and we could then detect very weak Mg n absorption 
doublets down to w r (2796) = 65 mA. 

We do not confirm the large strength of the z a b s = 1-0283 
Mgn syst em in PKS 0426-380 from the sample of Stocke 
& Rector i 19971 1. The spectrum of this Blazar is presented in 
Appendix C. 

In one Blazar, SDSS J094257.8-004705, there is an unusual 
broad absorption line with a FWHM 6000 km s _1 that we 
identify as a highly detached Mg n BAL at z a bs = 0.929 or Av = 
0.20c. More information on this system is given in Appendix 
B. Known and potential damped Lya systems are discussed in 
Appendix D. 

3.3. Statistical sample definition 

Our Blazar sample comprises 45 targets. Three are northern ob- 
jects without FORS1 data, but included here as they are in the 
Uy sample discussed by Stocke & Rector ( 11997l >. The S/N of 
their spectra only enables the detection of strong (w,(2796) > 
1.0 A) Mgn doublets. For the Blazars with unknown z em , the 
lower limit on z em is set by the Mgn system of highest z a b S . 
Associated Mgn absorbers, within 5000 km s of z em (or its 
lower limit) are excluded from our Mg n samples. 

The redshift path for all the sightlines is given in Table [3] 
and shown in Fig. [T] For both grisms, the minimum values Zmin 
correspond to either the wavelengths at which the overall effi- 
ciency is half that at peak efficiency or, if larger, the redshift 
at which the S/N enables the 3<x detection of the Mgn doublet 
(with w r (2796) limits of either 0.3 or 1.0 A). For the former, 
these values equal z min = 0.350 and 0.623 for the B600 and V600 
grisms, respectively. The efficiency of the B600 grism strongly 
declines towards its blue end and, although one weak Mg n sys- 
tem at z a b s = 0.2527 was detected in one sightline, the statistics 
are too poor to include this z range in our study. Similarly, the 
S/N decreases at the red end of the V600 grism (see Fig. A.l) 
due to strong sky lines which prevents detection of weak Mg n 
system at z a b s > 1.430 for most sightlines. This value is adopted 
as the upper limit of our z range for weak Mg n systems (the 
corresponding limit for strong systems is z a bs = 1.579). Thus 
the weak Mg n system at z a b s = 1 .5070 and the strong system at 
Zabs = 1-6832 (special setting: R600 grism) given in Table|2]are 
excluded from our statistical samples. 

4. Incidence of Mg n absorbers in Blazar spectra 

4. 1 . Global excess 

We measure the incidence of Mg n absorbers using our statistical 
sample of 45 Blazar sightlines. The number density per unit red- 
shift is dN/dz = N bs/Az, where N b s is the number of observed 
absorption systems within the redshift path Az. 

We first consider the entire Mgn absorber sample, i.e. all 
systems with w r (2796)>0.3 A. One of the strong system given in 
Table|2]is excluded (northern Blazar) since the low S/N prevents 
detection of weak systems. We thus get: 



dN . 
~dz 



(w r (2796) > 0.3 A) 



31 
25.11 



1.23!°;|at<z) 



0.83. (1) 



The errors are based on Poisson statistics for small numbers with 
limits corresponding to lcr confidence level of a Gaussian distri- 
bution as tabulated by Gehrels ( 1 19861 1 and also adopted by Tejos 
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Fig. 2. Left: Cumulative number of w r (2796)>0.3 A intervening Mgn systems towards Blazars, excluding the systems with Av < 
5000 km s _1 (thick, black curve). The curve of intermediate thickness (red) shows the cumulative number of such Mg n systems 
towards QSOs, adopting the incidence of these systems given by Nestor (see text), and the thinner (green) one is obtained by 
normalizing the latter to the total number of Mgn systems towards Blazars. Center: incidence of strong (w r (2796)>1.0 A) Mgn 
systems adopting the incidence of these systems given by Prochter et al. d2006a|). Right: incidence of weak (0.3 < w r (2796) < 1 .0 A) 
Mgn systems. 



et al. (12.009b . To compare this value to the incidence of Mg n ab- 
sorbers measured along QSO sight lines, we use the incidence 
of these systems per unit redshift derived from SDSS data. For 
w r (2796) > 0.3 A the incidence measured in QSO spectra is 
given by Nestor from his analysis of SDSS DR4 (private com- 
munication): 



dN 
~dz 



: 0.288 + 0.4869z- 0.091V 



(2) 



In Fig. |2] (left panel), we present the cumulative number of 
Mgn Blazar systems versus redshift for the entire sample. We 
also plot the cumulative distribution expected for Mg n QSO sys- 
tems (computed using Eq. 2) and then adjust the normalization 
for comparison with Blazars. There is no evidence for a differ- 
ence in functional form between Blazars and QSOs. 

Assuming that the Mg n Blazar absorbers are distributed at 
random along sightlines, as demonstrated for Mgn QSO ab- 
so rbers w ith Av > 5000 km s" 1 (Steidel & Sargent [1992l Wild et 
al. 120081 ). the probability distribution for N, the number of sys- 
tems expected in our set of sightlines, should follow a Poisson 
distribution: 



P(N obs > A^p) = Yj 

N=N 



N! : 



(3) 



where fi is the mean cumulative number expected for QSO sys- 
tems over the whole set of Blazar sightlines (the value is given 
by the endpoint of the red curve in Fig. |2j. With fi = 15.98 and 
No = 31, we find P = 0.00055, indicating that similar dN/dz for 
Blazars and QSOs are ruled out at a confidence level of 99.9%. 

The excess of Mg n Blazar systems with w, (2796) > 0.3 A 
equals: 



Ex -( dN ) /( dN ) 

V* /Blazar I * /QSO 



= 2.0+ 



0.4 
0.3' 



(4) 



therefore providing a new piece of evidence that the observed 
incidence of Mg n absorbers depends on the type of background 
source. This value is similar to that reported for GRB lines-of- 
sight. 



4.2. Dependence on rest equivalent width 

In order to investigate the origin of the excess, we measure the 
Mg ii incidence towards Blazars for strong (w, (2796) > 1.0 A) 
Mg ii systems. We get: 



dN 13 



0.43 ±° \ 6 2 at(z) = 0.82 . 



dz 29.93 u.iz— ( 5 > 

We adopt the fit given by Prochter et al. (|2006a) from their 
analysis of SDSS DR4: 

— = -0.026 + 0.374z - 0.145r + 0.026z 3 . (6) 
dz 

While we detect N Q b s = 13 systems, the number expected for 
QSO sightlines with the same z coverage is only 5.88. Using 
Eq.|3]with fx = 5.88 and N = 13, we get P = 0.0075, indicat- 
ing that similar dN/dz for Blazars and QSOs are ruled out at a 
confidence level of 99%. 

The excess of strong Mg n Blazar systems equals: 



E*(w r (2796) > 1.0 A) = 2.2 



,0.8 
fc 0.6 



(7) 



This excess is similar to that found for GRBs (Tejos et al. 2009; 
Vergani et al. 2009). In Fig. |2](center panel), we show the cumu- 
lative redshift distribution for the strong Mg n Blazar systems. 
We also plot the distribution expected for Mg n QSO systems 
(Prochter et al. 2006a) and that obtained after normalization to 
the total number of systems observed towards Blazars. 

For weak (0.3 < w r (2796) < 1.0 A) Mg n systems, we get: 



0.76 



t° f 7 at <z> = 0.83 . 



(8) 



dN _ 19 
~dz ~ 25.11 

Using Eq. [3] with No = 19 and the number expected for QSO 
sightlines with the same Az, fi = 10.83, we get P = 0.015. Thus 
similar dN/dz for Blazars and QSOs are ruled out at a confi- 
dence level of 98.5%, lower than for strong Mgn systems (the 
incidence of weak Mg n QSO systems is given by the difference 
in the values of dN/dz at (z) = 0.83 obtained with Eqs. 2 and 6). 
The excess of weak Mg n Blazar systems then equals: 



Ex(0.3 < w r (2796) < 1.0 A) = 1.7 ±£ 



0.5 



(9) 
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Table 3. Redshift path for strong (w r (2796) > 1.0 A) and weak 
Mgn (w r (2796) > 0.3 A) absorbers and S/N per pixel' 1 around 
peak efficiency. 



target 




z path 
Strong 


z path 
weak 


S/N 


0100-3337 


0.875 


0.494 


0.306 


50 


0210-5101 


1.003 


0.620 


0.620 


95 


0217+0144 


1.715 


0.956 


0.807 


190 


0234-3015 


1.690 


0.956 


0.807 


120 


0238+1636 


0.940 


0.558 


0.558 


80 


0241+0043 


0.989 


0.607 


0.607 


55 


0259+0747 


0.893 


0.511 


0.145 


40 


0334-4008 


1.351 


0.689 


0.689 


95 


0423-0120 


0.915 


0.533 


0.533 


180 


0428-3756 


1.110 


0.725 


0.725 


270 


0457-2324 


1.003 


0.620 


0.620 


350 


0538-4405 


0.890 


0.514 


0.514 


250 


0745-0044 


0.994 


0.611 


0.611 


200 


0811+0146 


1.145 


0.759 


0.759 


85 


0826-2230 


>0.910 


0.528 


0.528 


200 


0909+0121 


1.022 


0.634 


0.634 


140 


0942-0047 


1.362 


0.698 


0.698 


105 


0948+0839 


1.489 


0.812 


0.812 


55 


1009+0145 


> 1.085 


0.427 


0.191 


40 


1011+0106 


1.275 


0.614 


0.614 


60 


1132+0034 


1.234 


0.563 


0.563 


130 


1147-3812 


1.049 


0.665 


0.665 


55 


1245+0228 


> 1.096 


0.711 


0.703 


75 


1252-3319 


0.856 


0.475 


0.000 


28 


1338+1153 


>1.589 


0.923 


0.807 


150 


1408-0752 


1.494 


0.829 


0.829 


100 


1410+0203 


1.256 


0.595 


0.595 


75 


1419+0445 


> 1.684 


0.956 


0.807 


80 


1427-4206 


1.522 


0.857 


0.807 


100 


1522-2730 


1.294 


0.633 


0.633 


110 


1545+0032 


1.051 


0.394 


0.394 


50 


1743-0350 


1.054 


0.670 


0.670 


70 


1956-3225 


1.242 


0.582 


0.582 


50 


2015-0137 


>0.940 


0.558 


0.558 


1 15 


lAI 1 i 111(1 

2031+1219 


1.215 


0.555 


0.498 


55 


2134-0153 


1.285 


0.624 


0.624 


130 


2136-4443 


>0.980 


0.597 


0.597 


115 


2207+0041 


1.892 


0.817 


0.467 


55 


2214-2932 


1.636 


0.956 


0.807 


48 


2225-1113 


0.997 


0.614 


0.176 


35 


2225-0457 


1.404 


0.741 


0.741 


125 


2311-1038 


1.529 


0.864 


0.807 


75 


0221+3556 


0.944 


0.482 


0.000 


N/A 


0508+8432 


>1.340 


0.871 


0.000 


N/A 


1310+3220 


0.997 


0.534 


0.000 


N/A 


8 FWHM = 


3.4 px. 









Finally we note that our analysis indicates a possible differ- 
ence in the redshift evolution of the excess of Blazar and GRB 
strong absorbers. The excess of strong systems is higher at high 
Zabs along Blazar sightlines (a linear fit yields dN/dz ~ 1 .Oz-0.4) 
whereas the opposite trend is seen towards GRBs (Vergani et al. 
120091 1. No difference is seen for weak systems. The evolutions in 
redshift are shown in Fig.[3]for strong and weak Mg n systems. 



73 

N 
73 

\ 
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0.1 



N 
\ 

Z 
T3 



0.1 =- 
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0.5 



1.5 



Fig. 3. Number density evolution of strong, w r (2796) > 1 .0 A 
(top panel), and weak, 0.3 < w,(2796) < 1.0 A (bottom panel), 
Mg ii systems towards Blazars (blue squares), GRBs (red small 
triangles) as reported by Vergani et al. (120091 1 for their overall 
sample (with a downwards correction for the absolute values 
of strong systems: Vergani, private communication) and SDSS 
QSOs (solid line) as given by Prochter et al. (2006) for strong 
systems and Nestor (private communication) for weaker sys- 
tems. 

4.3. Results for BL Lacs only 

Since 35 out of the 45 targets in the overall Blazar sample belong 
to the BL Lac class, we can determine the incidence of Mg n BL 
Lac systems only, with still statistically meaningfull results. 

As in Sect. 4.1, we first consider the entire Mgn absorber 
sample (24 systems with w r (2796) > 0.3 A). This yields: 

24 

^-0.26 



^(w r (2796)>0.3A) = 
dz 18.33 



= 1.31^at(z)-0.84. (10) 



For the 10 strong (w r (2796) > 1.0 A) and 15 weak (0.3 < 
w r (2796) < 1.0 A) Mg n BL Lac systems separately, we get: 

10 

^0.14 



dN , 

22.81 



= 0.44 ±° n ] 9 A at(z) = 0.82 . 



(11) 



dN. 
~dz~ 



(weak) 



15 
18.33 



0.82 ±°ii at (z) = 0.84 . 



(12) 



These results are nearly identical to those obtained for the whole 
Blazar samples (Eqs. 1, 5 and 8), although with larger uncertain- 
ties. We thus find no evidence for a distinct behavior of BL Lac 
and FSRQ sightlines, which a posteriori justifies considering the 
Blazar sample as a whole. 

The excess of Mg n BL Lac systems equals: 



Ex = 2.1+^,2.2+ 



0.9 
t 0.7' 



1. 



,0.6 



(13) 



for the entire, strong and weak samples, respectively, while 
Stocke & Rector d!9971 > found on the basis of ten sightlines an 
excess of about 4-5 (note that among their five strong systems 
one is associated and another one is not confirmed). 
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4.4. Relative velocity distribution 

We now test whether the excess of strong Mg n Blazar systems 
arises from preferential values of relative velocity between the 
Blazar and the absorber. We consider the cumulative distribution 
of Mg n systems as a function of f3, 

v = (1 + z em ) 2 -(1 +Zabs) 2 (14) 

~C (1 +Z em ) 2 + (1 +Zabs) 2 ' 

where v is due to the expansion of the Universe and the peculiar 
velocity of gas ejected by a source at z em . 

In Fig. [4] we give the cumulative /3 distribution for strong and 
weak systems respectively. For comparison, we need to compute 
the distribution expected in /?-space from the background popu- 
lation of intervening Mg n absorbers. To this aim, we consider 
the redshift distribution adopted in Sect. 4. 1 and build the f3 dis- 
tribution for each sightline (due to varying z em values from one 
target to another, the transformation z — » f3 is specific to each 
sightline). Next, for any [3 bin, we combine all sightlines con- 
tributing to this bin. As expected, the same excess of systems as 
in Fig. |2]is apparent. More significantly, after scaling to the same 
total number of systems to assess whether the shapes of the dis- 
tributions are similar, one can see that strong systems display a 
relative excess of low /3 ~ 0.06 to 0.18 systems while the shapes 
of the QSO and Blazar distribution for weak systems are in quite 
good agreement. To estimate the significance of the excess of 
strong systems at low f3, we perform a KS test and find that we 
can reject the null hypothesis of a similar underlying distribu- 
tion only at the 80% significance level. This result must then be 
considered only as a trend. 

4.5. Associated systems 

Our results allow us to compare the number of Mg n systems 
at Zabs ~ z em to that found for QSOs; the latter is now well 
constrained by studies based on SDDS data (Wild et al. I2008t . 
Associated systems in Blazars might be more numerous than in 
QSOs due e.g. to a richer environment or to differing jet proper- 
ties. 

At least two associated systems with w r (2796) > 0.3 A 
are present; this is a lower limit because of the sightlines for 
which no emission redshift is available. The number density of 
w r >0.3 A QSO systems at z = 1 is dN/dz = 0.68 which implies 
an expected number of associated systems per sightline of 

dN dz 

AN= — I — | A/3 = 0.023, (15) 
dz d[3 

using 

|{||-(l+z) = 2. (16) 
dp 

Thus, for the 41 sightlines in which we can detect [3 - systems, 
we expect ~0.9 system. Towards radio-quiet QSOs, the fraction 
of w r (2796) > 0.5 A systems at Av < 5000 km s _1 , relative 
to the background level, is 1.7 while this factor reaches 3.3 to- 
wards radio-loud QSOs (Wild et al. 120081 Wild, private commu- 
nication: updated results based on the SDSS DR5 that are fairly 
insensitive to the w r limit). Thus, if Blazars have a similar excess 
of associated systems than radio-loud QSOs, we expect 3.3 asso- 
ciated systems at Av < 5000 km s~ l in our sample whereas there 
are two associated and three possibly associated systems. The 
latter could then be truly associated systems that indeed define 
the BL Lac/Blazar emission redshifts. 



5. Discussion 

5. 1 . Comparison between Blazars, GRBs and QSOs 

Our analysis provides a new piece of evidence that the observed 
incidence of Mg n absorbers depends on the type of background 
source. Towards both Blazars and GRBs, z ~ 1 Mg n absorbers 
are found in excess by roughly a factor of 2 with respect to QSO 
sightlines. In each case, the difference is significant at the ~ 3cr 
level. 

The reported excesses for GRB and Blazar sightlines are 
similar for Mg n absorbers with Wo > 1 A. For weak systems, 
comparison of the results obtained for Blazars and GRBs is 
less straightforward due to small number statistics for GRBs. 
Whereas we find an excess towards Blazars equal to 1.8±q5 at 
(z) = 0.83 (2.5<r significance level), there is no significant ex- 
cess for GRB absorbers (Vergani et al. 120091 1 with a value of 
1.2+Q4 (using Poisson statistics for the errors) at (z) = 1.57. 
Nevertheless, the GRB results differ from the Blazar ones at only 
a lcr significance level and it can be seen in Fig. [3] that the GRB 
absorber excess at their highest z bin follows the trend found at 
lower z for weak Mg n Blazar absorbers. It is thus not really pos- 
sible to conclude whether or not there is a difference between the 
weak Mg n Blazar and GBR absorber populations. 

Our analysis also suggests a possible difference in the red- 
shift evolution of this excess but larger samples are needed to 
make definite statements. 

5.2. Dust and amplification biases 

Dust extinction and/or gravitational magnification due to fore- 
ground matter can modify the observed density of background 
sources and therefore the observed incidence of intervening ab- 
sorbers. Several authors already pointed out these effects as po- 
tentially responsible for (some of) the Mg n absorber excess in 
GRB spectra (Porciani et al. 2007). Here we explain why these 
effects are unlikely to have a significant contribution in giving 
rise to an apparent excess of strong Mg n absorbers both in GRB 
and Blazars spectra compared to the value measured with QSOs. 

GRBs are preselected on the basis of their gamma-ray flux. 
Rapid follow-up observations of the afterglow emission allow 
us to observe these objects while their apparent magnitude is 
significantly (i.e. several magnitudes) brighter than typical lim- 
iting magnitudes. This difference is much larger than the typical 
levels of extinction induced by intervening Mg n absorbers. The 
dust content of these systems has been quantified by Menard et 
al. J2008 ): these authors showed that on average, Mg n absorbers 
with w,-(2796) ~ 1 A give rise to a color excess E(B - V) - 0.01 
mag. These authors estimated that about one percent of SDSS 
QSOs are missed due to the presence of dust in foreground Mg n 
absorbers. Most BL Lac objects are identified on the basis of 
their radio or X-ray emission (e.g. Padovani & Giommi |l995l ) 
or their gamma-ray emission. This preselection is not sensitive 
to a dust obscuration bias (only low frequency X-rays would be 
affected by absorption but very large hydrogen column densities 
are required). The optical follow-up could be affected by such 
an effect. However, Blazars are variable in this wavelength range 
with a typical amplitude reaching about one magnitude, and for 
a few Blazars of our sample up to 4-5 magnitudes. Consequently, 
as in GRBs detection of Mg n absorbers towards Blazars is not 
affected by dust extinction. 

Gravitational magnification can change the density of back- 
ground sources. This effect depends on the gradient of the mat- 
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Fig. 4. Left: Cumulative number of strong (w r (2796) > 1 .0 A) intervening Mg n systems towards Blazars versus fi = v/c, excluding 
the systems with Av(z e m,Zabs) < 5000 km s _1 or /3 < 0.0167 (thicker, black curve). The curve of intermediate thickness (red) 
shows the predicted number of strong Mg n systems, adopting the incidence of these systems towards QSOs derived by Prochter et 
al. d2006al >. and the thinner (green) one is obtained by normalizing the latter to the total number of strong Mg n systems towards 
Blazars. Right: Same for weak (0.3 < w r (2796) < 1.0 A) Mgn systems. 



ter over-density along the sightlines and the slope of the source 
number counts as a function of magnitude (Naravan ll989l >: 

n(m) = ff™- 1 n (m) , (17) 

where n and no are the lensed and unlensed density of sources in 
the sky, /i is the gravitational magnification along their sightline 
and a(m) is the slope of the source number counts as a func- 
tion of magnitude. The average magnification effect induced by 
galaxies was detected by Scranton et al. (f2005). These authors 
report an average magnification excess of order one percent for 
impact parameters of about 50 kpc. Gravitational lensing by ab- 
sorbers has been theoretically addressed by a number of authors 
(Bartelmann & Loeb |1996l Menard EUUBl Given the typical im- 
pact parameters of strong Mgn absorbers (~ 20-50 kpc) the 
magnification excess fi - 1 is expected to be small, i.e. of or- 
der a few percents. Observationally, constraints on QSO magni- 
fication by Mg ii absorbers have been reported by Menard et al. 
(2008). These authors report /i - 1 < 0.1. This method cannot 
be applied to the Blazar population since their number counts as 
a function of magnitude is ill defined due to their strong opti- 
cal variability. But, the latter being far larger than any expected 
gravitational magnification, any significant amplification bias is 
most unlikely for Blazars. Finally, HST imaging surveys show 
that gravitational lensing appears to be unimportant to the BL 
Lac phenomenon (Giovannini et al 2004|. Recently Wyithe et al. 
( 120101 1 proposed strong lensing as an explanation for the higher 
incidence of Mg n absorbers in front of GRBs. They argue that 
20 to 60% of the GRBs with strong Mg n absorbers are strongly 
lensed. Here we note that if this was the case, the GRB host 
galaxy should also be strongly lensed and appear as arcs or mul- 
tiple images. So far HST imaging of GRB host galaxies have not 
reported such a phenomenon (Conselice et al. 12005 l l. indicating 
that strong lensing is not often induced by intervening Mg n ab- 
sorbers. 

Therefore, while dust extinction can lower the apparent inci- 
dence of absorbers in front of quasars and gravitational lensing 
can increase it in front of GRBs and Blazars, the amplitude of 
these effects is at least an order of magnitude lower than what 
would be required to explain a factor ~2 change in the incidence 
of Mg ii absorbers. 



5.3. An intrinsic origin? 

We now consider whether a significant fraction of the Mg n sys- 
tem excess could be intrinsic by exploring some of the physical 
conditions required for cold gas to be entrained at high veloci- 
ties. Such a scenario is appealing in the context of Blazars since 
these sources are known to display relativistic jets pointing to- 
wards us. If material is present in the neighborhood (either gas 
surrounding the AGN, within the interstellar medium of the host 
galaxy, or further away within its halo) it might be swept up by 
the jet and accelerated to high velocities. In our sample, there 
is no marked accumulation of systems at /3 < 0.02 thus, if such 
systems give a significant contribution (^ 50%) to the observed 
Mg ii systems, ejection velocities of a few x 10 4 km s must be 
involved. 

Towards QSOs, outflow velocities in this range are currently 
observed for BAL systems. The latter are thought to arise when 
the line of sight intersects the wind emanating from the accre- 
tion disk with the appropriate inclination angle - typically 30 deg 
(see Elvis 2000) - which accounts for the large velocity disper- 
sion found in these systems. Regarding Blazars, the inclination 
is supposed to be smaller (around 10 deg; Fan et al. 120091 ) and as 
a result, the line of sight is not expected to cross the disk wind. 
Indeed, BAL absorption is usually not seen in Blazar spectra. 

One may rather invoke material impacted by the jet along its 
way as the source of the absorption. This could be for instance 
diffuse gas located relatively close to the AGN. In such a config- 
uration, jet material is expected to sweep up the external gas and 
while the corresponding column density (N e ) increases, the ad- 
vance velocity of the interface (V) decreases (Falle [19911 Kaiser 
& Alexander [19971 Komissarov & Falle [199"8l ). If we identify 
the advance velocity with the apparent ejection velocity of the 
absorbing gas, the question is then: can we still have total col- 
umn densities of swept up gas as large as 10 18-20 cirT 2 together 
with velocities as high as a few 10 4 km s _1 ? Komissarov & Falle 
( 119981 ) have shown that in the early phase of a relativistic jet 
growth V c while later, the dynamics of the interface follows a 
power-law behavior with respect to the time t, with the jet length 
(L) increasing as f 3 ' 5 , corresponding to a velocity V = dL/dt 
decreasing as r 2 ^ 5 . In this regime, one can use the power-law 
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expressions provided by Komissarov & Falle d 19981) to relate di- 
rectly L to V (by eliminating f), which leads to 



(18) 



The column density of swept up gas, N e , simply writes N e = 
L x n e . Assuming for simplicity a pure hydrogen gas to connect 
p e and n e , we finally get the numerical relation: 



N e (cm- 2 ) 2xl0 19 



I0 45 ergs 



lO.lcm- 3 / VO.lc/ 



This relation indicates that for reasonable jet powers and ambi- 



ent gas densities, a column density as large as 10 u cm can 
be swept up, with a relative velocity attaining 0.1c. The above es- 
timate for N e is an upper limit for the amount of gas effectively 
present at the interface because part of the swept up external gas 
flows away from the contact zone. However, even if only a small 
fraction (e.g. 0.1) remains at the interface, the column density 
can be large enough to induce the observed absorption systems. 
Note that such an intrinsic scenario does not necessarily imply a 
large velocity dispersion since, contrary to the case of QSO BAL 
systems, the absorbing material is not closely associated with the 
active nucleus. 



6. Conclusions 

Attempting to shed light on the unexplained excess of Mg n ab- 
sorbers in front of GRBs, we measured the incidence of such 
absorbers in front of another class of objects: the Blazars. 

We observed with FORS1 at the ESO-VLT 42 Blazars with 
an emission redshift 0.8 < z em < 1.9, to which we added the 
three high z northern objects belonging to the Uy BL Lac sam- 
ple. From this sample we detect 32 Mg n absorbers in the redshift 
range 0.35-1.45. Our main results are as follows. 

1 . We find an excess in the incidence of Mg n absorbers com- 
pared to that measured towards QSOs by a factor ~ 2, de- 
tected at 3cr. The amplitude of the effect is similar to that 
found along GRB sightlines. 

2. In front of Blazars, the excess is apparent for both "strong", 
i.e. w r (2796) > 1.0 A, and weaker systems, with 0.3 < 
w r (2796) < 1.0 A. 

3. The dependence on velocity separation with respect to the 
background Blazars indicates, at the ~ 1 .5cr level, a potential 
excess for 0.06 </3< 0.18. 

This work brings another piece of evidence that the incidence 
of extragalactic absorbers systems might depend on the type of 
background source considered. We argue that biases involving 
dust extinction, weak and/or strong gravitational amplification 
are not likely to affect the frequency of Mg n systems detected 
towards either Blazars and GRBs at a significant level. We also 
briefly discuss the physical conditions required for these ab- 
sorber systems to be physically associated with the Blazars. 

We claim that a satisfactory explanation for the apparent ex- 
cess of Mg n absorbers in front of certain types of background 
sources still awaits. Further investigations should be made on 
both theoretical and observational sides. For the former, realistic 
numerical modelling of the jet/ambient gas interaction should fo- 
cus not only on jet properties as commonly done (see e.g. Choi 
et al. 12007b but rather on those of the swept up gas at differ- 
ent stages of the jet growth. We stress in particular that taking 
into account the thermal behavior of the gas (heating/cooling 



processes) is important to properly assess its physical properties 
(ionisation level, spatial structure, temperature: e.g. Mellema et 
al. l2002l l. 

Given the lack of current theoretical understanding, more 
observational results might offer us some guidance, especially 
by exploring dependencies with source properties. A first test is 
to split QSOs in various sub-classes. The most straightforward 
one involves their radio brightness: this has been investigated by 
Ellison et al. £2004) who have shown that the incidence of Mg n 
absorbers does not appear to depend on the radio loudness of 
the background QSO. Other tests concern the size of the opti- 
cal emission regions. For QSOs, Pontzen et al. (|20O7J reported a 
consistent Mg n incidence measured over the continuum regions 
and the larger, C m emission line regions. Such a test should also 
be conducted for Blazars in high and low states, when the contin- 
uum flux originates mainly either from the jet or from the accre- 
tion disk, respectively. Those Blazars from our sample display- 
ing several systems are obvious targets for such a study. Repeat 
high resolution spectra would in particular constrain possible 
variability of various absorption lines (Mg n, Mg i, Fe n), thus 
test the intrinsic origin scenario, and yield kinematics and ion- 
ization properties of Mg n Blazar absorbers to be compared with 
those of Mg ii QSO absorbers. 
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Appendix A: New emission redshifts 

The spectra of three Blazars show very weak emission lines yielding new emission redshifts. They are presented in Fig. I A. 1 
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Fig. A.l. Blazars with new emission redshift (Mg n and C n] emission lines). 



In PKS 0332-403 (Fig. lA.ll top panel), the weak emission at A — 6580 A is identified as Mgn at z em = 1.351. As reported in 
Table[2] there is a weak Mg n absorption doublet at z a bs = 1 -2083 with associated Fe n absorptions (the /12600 doublet and the /12382 
triplet); the other Mgn system at z&s = 1-0791 is too weak to be included in our statistical Mgn sample. 

In PKS 1008+013 (Fig. IA.ll central panel), there is an asymmetric emission line centred at A — 6368 A identified as Mgn at 
Zem = 1-275, which is also clearly present in the SDSS spectrum of this source. We thus do not confirm the 'uncertain' redshift, 
Zem = 0.8615, given by Collinge et al. ( 120051 . We do not detect any absorption line in the FORS1 spectrum of this Blazar. 

In PKS 1407+022 (Fig. IA.ll bottom panel), there are two weak emission lines at A — 5238 and 6308 A that we identified as 
C n] and Mg n at z em = 1 -253. There is a weak Mg n absorption doublet at z a bs = 1 ■ 1 123 (see Tabled with possible associated Fe n 
absorptions (marginal detections of Fe n2382 and Fe n2600). 



Appendix B: The discovery of a broad absorption line in a Blazar spectrum 

For one Blazar (the weak FIRST radio source SDSS J094257. 8-004705), there is an unusual broad absorption line with a FWHM 
- 6000 km s _I as shown in Fig. IB. II This feature cannot be the signature of a dust shallow absorption at 2175 A associated with 
the Blazar (z em = 1-362), since it should be at A r > 2285 A. Thus, it is identified as a Mg n BAL at z a t, s = 0.929 or Av = 0.20c. Such 
highly detached BALs, although very rare, have already been detected at low redshift (e.g. the high ionization BAL in PG 2302+029 
at Zabs = 0.695 reported by Jannuzi et al. d!996l >). In sources at 1.5 < z < 3 from the FIRST Bright Quasar Survey (FBQS), there is 
a 5% fraction of low ionization (Lo) BALs (Becker et al. 2000); the frequency of LoBAL QSOs is dependent of radio-loudness and 
decreases for the most radio-luminous QSOs. We thus expect at most two LoBAL QSOs in our Blazar sample since 33 of the 45 
Blazars are strong (~1 Jy) radio sources, which is consistent with our single LoBAL detection. 
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Fig. B.l. BL0942-0047 with C iv and Mg n emission lines and a Mg n BAL at A ~ 5400 A. 

There are also two narrow absorption systems in the Blazar spectrum. The strong Mg n doublet at z a b s = 0.8182 has Mg i, Mn n 
and Fe n associated strong absorptions. The weak system at z a b s = 1.0231 shows Mg n absorption only. 



Appendix C: Not confirmed strong Mg n absorption in one source of the Stocke & Rector BL Lac sample 



PKS 0426-380 is one of the ten sources presented by Stocke & Rector d 19971 1. sample for which an excess of strong Mg n in BL Lac 
objects had first been reported. The high S/N (270) of our FORS1 spectrum enables us to detect very weak absorption lines. The 
absorber at z a b s = 1.0283 is not confirmed as a strong Mg n system (w r (2796) = 0.56 A) and has weak Fe n(2382,2600) and Aim 
doublet associated absorptions. 

There is a stronger Mg n doublet (w r (2796) = 0.93 A) at z a b s = 0.5592 with associated Mg i absorption together with the five 
stronger Fen lines. We also detect Galactic Can absorption, as well as the only intervening (z a b s > 0) Can absorber in our Blazar 
sample. The latter is at z a b s = 0.1940 with a marginal detection of the weaker line of the doublet but exactly at the same redshift. 
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Fig. C.l. In this BL Lac the Mg n absorption at z a b s = 1 .0283 is weaker than reported by Stocke & Rector dl 9971 1 . 



Appendix D: Known and potential DLAs 

Classical Damped Lya- systems (DLAs) have H i column densities N(H i) > 2 x 10 20 cm 2 . The fraction of DLAs among intervening 
QSO absorbers, with constraints on Mg n2796 and Fe n2600 (or the ratio of their equivalent widths) and on Mg i2852 is about 
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40% (Rao et al. 120061) . Is the DLA fraction similar for Blazars? Among our sample of 45 Blazars, N(Hi) can be derived from Lya 
absorption line profiles for eight sources only, i.e. those with HST/FOS or STIS spectra. Since most of our selected Blazars are 
strong radio sources, the DLA subpopulation which traces a cold gas phase can also be identified by 21 cm absorption. 

Three of the strong Mg n Blazar absorbers have UV spectra. One is a DLA, the z abs = 0.5245 absorber towards PKS 0235+164, 
and two are sub-DLAs with N(Hi) < 1 x 10 20 cirT 2 , the intervening absorber at z abs = 1-3439 towards PKS 0215+01 and the 
possibly associated absorber at z abs = 0.9107 towards PKS 0823-223. The other properties of DLAs and sub-DLAs mentioned 
above are their large w r (Fe n2600)/w r (Mg n2796) ratio (>0.5) and Mgi2852 strength (e.g. Bergeron & Stasirika [19861 Rao et al. 
2006). Four strong Mgn absorbers do not satisfy these criteria and thus cannot be DLAs: (i) two intervening ones at z abs = 0.7753 
in SDSS J024156.4+004351 and z abs = 1.4247 in SDSS J094827.0+083940, (ii) two possibly associated systems at z abs = 1-0851 
in PKS 1008+013 and z abs = 1.6832 in SDSS J141927.4+044513. 

Two strong Mg n systems in the BL Lac sample of Stocke & Rector ( 1997), also in our sample, were known 21 cm intervening 
absorbers, at z abs = 0.5245 (DLA mentioned above) and z abs = 0.6850 towards the lensed source B2 0218+357. There are three 
other intervening, strong Mg n absorbers towards bright radio sources, with DLA signatures (detection of weak transitions of Fe n, 
Znn and Cm). We have recently obtained GMRT and GBT time to observe these sources. The GMRT data show 21 cm absorption 
atz abs = 1.2753 in PKS 1406-076 and z abs = 1.1158 in PKS 2029+121 (Gupta et al. in preparation). The FORS1 spectrum of PKS 
1406-076 is presented in Fig. ID. II Analysis of the GBT data for PKS 0454-234 is in progress to search for 21 cm absorption at 
z abs = 0.8922. 

In our Blazar sample, 10 strong Mg n systems satisfy the criteria of high N(H i): four are confirmed DLAs and one is a probable 
DLA, thus yielding a DLA fraction of 40% compatible with that determined for QSOs absorbers. Additional UV and/or 21 cm data 
are needed to ascertain whether the incidence of DLA in our strong Mg n-Fe n Blazar sample is even higher or not than for QSOs. 
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Fig. D.l. This Blazar (PKS 1406-076) is a variable FSRQ, with superluminal motions, and also a gamma-ray emitter (see Sect. 
12. II and Table Q~|i. It has strong emission lines (Mgn and Cm]) and two Mgn absorption systems: a strong one at z abs = 1.2735, 
and a very weak one at z abs = 1.2913 thus not included in our statistical Mg n sample. The former is a 21 cm absorber and a DLA 
since, in addition to Mg i and the five stronger Fe n lines, the very weak transitions of Fe n2249,2260, Zn n2026,2062blend and 
Crn2056,2062blend,2066 are detected. 



